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q ■ the exact diagonalization of algebraically determined matrixes based 
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on powers of Morse coordinates is made substantially more efficient 
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by choosing a properly adapted quantum-mechanical basis, specifi- 
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cally tuned to the molecular potential. A substantial improvement is 
achieved while still retaining the full advantage of the simplicity and 
numerical light-weightedness of an algebraic approach. In the scheme 
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we propose, the basis is parameterized by two quantities which can be 
adjusted to best suit the molecular potential through a simple mini- 
mization procedure. 

Keywords: vibrational spectra, algebraic method, Morse oscillator, quasi number 
state basis, basis optimization, anharmonic vibrations. 

Introduction 

In a previous workP an algebraic method for the computation of vibra- 
tional spectra of diatomic molecules was introduced. Although this is a 
1-dimensional (ID) problem, thus an in principle trivial task, the algebraic 
method shows substantial advantages over both the real-space grid solution 
of the Schrodinger equation and harmonic-oscillator-based techniques. These 
advantages are especially important for extensions to the multidimensional 
problem of polyatomic vibrations. 

The expansion of the molecular potential in powers of the Morse-potential 
related quantity v(x) = e~ a ( x ~ x °) — 1, namely 



allows an efficient and accurate approximation of a well-behaved molecular 
potential in the whole energy range, from the minimum region to the dissoci- 





k=2 



ation threshold, generally involving a moderate number iV max + 1 parameters 
a2, • • • , at and xq. Even potentials substantially distorted with respect to the 
Morse potential can be treated successfully. With the potential expressed in 
the form of Eq. ([T]) , the complete Hamiltonian 



(here \x is the reduced mass of the 2-body problem and x is the radial coor- 
dinate) can be represented on a quantum-mechanical basis of choice. 

The accuracy and efficiency of the direct diagonalization methods rely 
both on the accuracy of the potential approximation of Eq. ([T]) and on the 
properties of the selected basis. The basis had better be complete but also 
manageable, i.e. related to the algebraic properties of v(x), so that the evalua- 
tion of the matrix elements can be done rapidly and without approximations: 
this will be needed especially in view of extensions to polyatomic molecules. 

2 The Basis 

Previous researchP^ showed that the basis 
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+ V d (x) 
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4>n(y) 
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y'eSltf-^y), a > 0, n = 0, 1, 2,... 
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with 

y(x) = (2s + 1) e - a{x - xo) , (4) 
can be usefully employed in general diatomic contexts, with the special choice 

<t = 8-[s], (5) 

where [s] indicates the integer part of s, and with s related to the Morse 
term a 2 (v(x)) 2 in the potential expansion jT]), by 



V2/ia 2 1 

S = — -r (6) 



With the conditions (loTUj) the basis was named quasi number state basis 
(QNSB). 2 In the present work, we only assume a and xq in Eqs. (I3f4"l) are 
the same as in the potential expansion (JT]), and that a > and s > — |, but 
release all additional unnecessary conditions on o and s, for example those 
expressed by Eqs. ( 151161) . or the condition defined by Tennyson and Sutcliffe^ 
(TS): 

with s fixed by Eq. (jBJ). Equation (jHJ) thus defines a (s, cr)-parameterized 
family of bases, generalized QNSB (GQNSB), all sharing the following main 
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features: (i) the basis © is complete; (ii) the kinetic and potential operators 
can be written in terms of generalized ladder operator as specified below, so 
that (iii) the matrix elements of a vast class of relevant operators is com- 
putable easily and exactly by means of simple algebraic relations.^ 

Even though all infinite GQNSB's are substantially equivalent, regard- 
less of s and a, different bases characterized by different values of s and 
a show different performances when truncated to a finite number N s of 
states and applied to a given quantum mechanical problem specified by 
fj,, a, 0*2,0,3, ... , ajVmaK- Indeed, the purpose of the present work is to demon- 
strate that a properly chosen truncated GQNSB can improve the efficiency 
of the computation substantially, compared to earlier choices)^ 

3 Matrix elements 

We follow here the same approach^ derived from SUSY quantum mechan- 
j cs pi] ^ e introduce the generalized Morse ladder operators^! 



A(q) 



qI ~2 + M P * 
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parameterized by the real quantity These operators, with a suitable choice 
of q, act on the states (EJ) of the GQNSB as ladder operators: 



A(o + n) <j>. 



•n 



(9) 



A\a + n) <j>, 



n 




where 



C n = y/n(n + 2cr-l) . 



(10) 



According to Eqs. flHIIUl) . a links the parameterized basis © to the corre- 
sponding family of generalized ladder operators. Thus, each and every basis 
of the form of Eq. ([3]) can be managed algebraically in this formalism, for 
any given choice of a > 0. In practice, the eigenfunctions depend ex- 
plicitly on s, a and xq, beside a. We fix to the position of the minimum 
of the potential (pQ), as it would not provide a substantial advantage other- 
wise. Likewise, we select for a the same value as in the potential expansion, 
because otherwise all relevant matrix representations would be dense rather 
than sparseP With these constraints on xq and a, an arbitrary s can be 
usefully employed in the basis definition: for any s value, the momentum op- 
erator p x and the multiplication operator e - a ( x - x o) can De written in terms 
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of the ladder operators 



-a(x-xo) 



Px 



2ql - 



A\q)+A(q) 



ha 
~2i 



(2s + 1) 



;n) 
;i2) 



where also the A operators depend implicitly on the s parameter appearing 
in the definition (J4]) of y. On the GQNSB ([3]), the matrix elements of any 
physical operator expressed as a polynomial of e~ a ( x ~ x °^ and p x can be com- 
puted algebraically since Eqs. (I11|I12I) express them in terms of the ladder 
operators of the corresponding specialized basis. We derive here explicitly 
the algebraic form of the Morse Hamiltonian for general q and s. 
Using Eq. ( fl2l) . the kinetic operator K = becomes 



2^,2 



K 



h a 



\A 2 (q) + i t2 (g) - A{q)A\q) - A\q)A{q)\ . (13) 



By applying the commutation relations 



[i(g),iV)l = (? + QV -(A(q) + ti (q% 



(14) 



[A(q),A(q')} = [A\q),A\q')} = 0. 
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K reduces to 



j2 2 

K = --^{A 2 (q) + AV(q) - 2ql + A(q) + A\q) - 2A\q)A{q)] . (15) 

OTfl 



Likewise, powers of e a ^ x appearing in the potential-energy operator are 
obtained starting from Eq. ( TTTT) . For example, 



e -2a(*-*„) = ( 16 ) 
1 



-{4q 2 f - iq[A\q) + A(q)] + A 2 (q) + A^(q) + A(q)A\q) + A(q)A(q)} 

(zs + ly 

(2g | i)2 {2(2g 2 + q)I - (Aq + l)[A\q) + A(q)] + 2A\q)A{q) + A 2 (q) + A^ 2 (q)} 



Thus, the Morse-potential term reads 



(^)) 2 = 77 r ^ T ^{2(2g 2 + g)/-(4g + l)[it( g ) + i( g )] + 2it( g )i( g ) 

(2s + ly 

+ A 2 ( q ) + A*{q)} - j^^iW - [Mq) + Mi)\} 

' -{(4g 2 -2q- 8sq)I + (As - Aq + l)[i f (g) + 1(g)] 



(2s + l) 2 

+ 2A\q)A(q)+A 2 (q)+AV(q)}. (17) 
Accordingly, the Morse Hamiltonian Hm = K + <X2(v(x)) 2 is expressed in 
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algebraic form as 



H 



M 



2a 2 



2^,21 



+ 



h a 
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a 2 



A\q)A(q) + q 



2^,2 



2a 2 . h a 

' (2q - 1 - 4s) + 



(2s + 1 



4m 



(2s+l) : 



(4s -4q + l) 



h 2 a 2 



5m 



[i+(g)+A(g)] 



a 2 



h 2 a 2 



(2s + l) 2 8m 



[A 2 (q) + AV(q)}. 



The representation of Eq. (1181) shows that the Morse Hamiltonian is generally 
5-band diagonal on a GQNSB of the form (jHJ). We stress that the expression 



(jTBj) holds for any choice of parameters s and g, regardless of them being 
connected to any specific physical constraint. 



If the condition 



d-2 



2^,2 



fi a 



(2s + 1) 



5m 



(19) 



(equivalent to Eq. (jSJ)) is satisfied, then the last term, proportional to [A 2 (q) + 
A^ 2 (q)] drops from Hm- In other words, the choice of the parameter s 
of Eq. ^ makes the Morse Hamiltonian tridiagonal on the corresponding 
GQNSB basis, irrespective of q. Under this special condition ( fl"9l) . the Morse 
Hamiltonian simplifies to: 



H 



M 



j^y 2 { [A\q) + A{q)\ (s - g) + A\q)A{q) + (q 2 - 2gs)/} . (20) 
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The form of Eq. f[2"Ul) . indicates that by further setting 

q = s, (21) 

the operator form of the Hamiltonian simplifies even more, and the Morse 
Hamiltonian factorizes as: 

Hm = 4 (^T^[ it (*) i ( s ) - ^ > ( 22 ) 

which recovers the algebraic form of the Morse Hamiltonian of previous 
works 

The use of different values of q and s produces a GQNSB, where the alge- 
braic computation of the matrix elements of the Hamiltonian (j2J) is not sig- 
nificantly more intricate: in particular on a GQNSB, the Morse Hamiltonian 
is 5-band diagonal, rather than tridiagonalP and higher powers of (v(x)) k in 
Eq. (0Q) generate (2k + l)-band diagonal matrices (like in the QNSB). 

For practical potentials, usually substantially distorted from the pure- 
Morse (v (x)) 2 term, the actual eigenfunctions can be represented poorly by 
the [s] + 1 Morse bound states, or equivalently by their QNSB counterparts: 
to achieve a good convergency of all eigenfunctions, the QNSB often needs to 
be complemented by a large number of states, far beyond [s] + 1. A suitably 
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chosen GQNSB can thus prove significantly more efficient, especially in a 
multi-oscillator polyatomic context. 

4 GQNSB parametric dependency 

The shape of the wavefunctions ([3]) depends on the four parameters 
Xq, a, s and a: different shapes imply different convergence properties when 
employed to build the matrix representation of the Hamiltonian. A brief 
analysis of the dependency of the shape of GQNSB states on the various pa- 
rameters can be useful to gain some insight in their role. Figure [T] shows the 
profile of three states of the form (jHJ), under conditions (jSJ) and (EJ). Note that 
the n = state is located substantially at the right of the Morse equilibrium 
position Xq, and that further states move in toward Xq for increasing n. This 
contrasts with the behavior of a basis of energy eigenstates of a well centered 
in Xq. Figure [2] illustrates the behavior of a GQNSB wavefunction, Eq. 
after variation of the parameters s and a involved relative to the QNSB val- 
ues, Eqs. (jSl EJ). The dependence on the s parameter (Fig. [2^) is weak: by 
increasing s, the eigenfunction shifts almost rigidly towards the outer region. 
The cr-dependence (Fig. [2b) is less trivial: for larger a, the wavefunction 
deforms and shrinks, concentrating toward the region of the minimum, and 
decaying more rapidly at large x. The role of the a parameter is particularly 
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important: as the n th GQNSB wavefunction has the general form 



(j) n (y) oc e 



*f 2 y°Pol\y,n]> 



(23) 



(Pol[y, n) stands for a polynomial of degree n in the variable y), a controls 
the decay rate of the wavefunctions for y — > 0, i.e. at the dissociation region. 
In particular, by choosing small a, the basis wavefunctions spread away from 
the well region thus improving the convergency of high-energy states, possibly 
at the expense of quality of the low-energy states in the well. Equation f l23l 
and Fig. [2] show that the general shape and in particular the amount of 
localization of the GQNSB wavefunctions can be tuned freely by choosing 
suitable s and a parameters: this allows improving the variational efficiency 
of a truncated GQNSB for a specific quantum-mechanical problem. 

5 Optimization of the basis parameters 

Assume that the exact iVj, bound state eigenvalues of the Hamiltonian 
are known; we can measure the RMS discrepancy of the discrete spectrum 
due to basis-incompleteness by 

, N b -1 




n=0 



(24) 
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in terms of the numerical eigenvalues E i: obtained by diagonalizing the ma- 
trix of H, Eq. (EJ), on a finite GQNSB composed by the first N s (> Nb) states 
and parameterized by s and a. For fixed N s we can search for the optimal 
s min and <r min that make A minimum. 

In fact, the a priori knowledge of the exact eigenvalues Ef x is not nec- 
essary: due to the variational nature of basis truncation, a "better" basis 
makes all eigenvalues E^ lower. Accordingly, the optimal s min and <r min pa- 
rameters can be defined as those producing the lowest eigenvalue spectrum 
for the assigned basis size N s , i.e. those minimizing 

, N b -1 

A = wX E - (25) 

6 n=0 

This approach only requires that the number Nb of bound eigenstates is 
known. Of course, the number Nb of bound states can be determined once 
and for all, for example by means of a calculation on a very extended QNSB. 
The minimization of A and of A leads generally to slightly different results, 
but the following qualitative discussion applies equally well to both schemes. 
Unless specified, for the determination of s min and o" min , we minimize A as 
defined in Eq. (|25l) . and compare A to its fully-converged value A computed 
on a largely complete basis. 

In a typical application of the GQNSB, one starts from a molecular po- 
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tential energy expressed in terms of an expansion of the form of Eq. ([T]). 
Before considering realistic dimers (H2 and Ar2), we illustrate the properties 
of the optimized GQNSB for a simple toy potential defined by 

N max = 4, a2 = 04 = 625, 03 = 0, a = 4 , and xq = 1, (26) 

which we solve combined with a kinetic term specified by % = 1, /i = 1. 
We minimize A with respect to s and a, for two fixed numbers of basis 
states N s = 30 and 16. Figure [31(a) shows the values of the individual eigen- 
value discrepancy (E n — -E° x )/a 2 for the potential f[2"6"]) . for the QNSB, for 
an optimized GQNSB (OGQNSB), and for s and a chosen according to the 
prescription of TSP^ The optimized parameters of the N s = 30 OGQNSB 
are s min = 20.01 and cr min = 0.435, to be compared with the QNSB ones 
s = 8.338 and a = 0.338, and those chosen according to the prescription of 
TS 31151 s = 8.338 and a = 9. For this potential A = -444.90, and the corre- 
sponding A — A are 3 ■ 10~ 6 for the OGQNSB (A equaling A to 5 decimal 
digits), 0.061 for the QNSB, and 457 for the TS choice. Both QNSB and the 
OGQNSB retrieve all the bound states, but the OGQNSB produces much 
better converged eigenenergies, especially near dissociation. The TS basis 
instead yields only 9 of the 14 bound states, only few of which are converged 
within 10~ 2 a2, which explains the large discrepancy A — A . 
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Figure [31(b) shows the same individual discrepancies obtained with a basis 
of N s = 16 states instead of 30. The s and a values of the QNSB and the 
TS basis are of course unchanged, while for the OGQNSB they change to 
■Smin = 14.47 and cr min = 0.314. The discrepancies A— A deteriorate to 0.106, 
87.35, and 4045.9 for OGQNSB, QNSB and TS respectively. Clearly the 
OGQNSB maintains a fair accuracy throughout the spectrum, by allowing 
for slightly less accurate lowest bound states, at the benefit of those near 
dissociation. In contrast, the N s = 16 QNSB fails in obtaining the two bound 
states closest to dissociation, and the TS basis only produces 6 bound states. 
Thus, basis parameters optimization allows a substantial improvement of the 
accuracy of the results, with the same computational cost. In other words, 
the convergence speed of the computation can be improved drastically by 
means of a suitable choice of s and a, for example Fig. [3] demonstrates an 
equal accuracy of the OGQNSB of 16 states and the QNSB of 30 states. 

Figure H] illustrates a typical s dependence of the total discrepancy A— A : 
for cr equal to its optimal value a m i n (solid curve), as s approaches the optimal 
■Vm value from below, A decreases relatively slowly, while for s increasing 
beyond s m j n , A grows very steeply. The a dependence of A has a sharp and 
roughly symmetrical deep minimum around cr min . 

The reason for the observed s and a dependencies of A is related to 
the GQNSB wavefunction profiles of Figs. Q] and [2} and Eq. (1231) . When s 
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increases, the GQNSB wavefunctions shift almost rigidly toward the disso- 
ciation region of the potential. As the GQNSB wavefunctions decay much 
more rapidly for small x than for large x, approaching s m i n from below the 
accuracy of the representation of the bound states localized in the well re- 
gion improves slowly, but soon after the optimal s is found, all wavefunctions 
move their localization region to the right of the equilibrium position, and 
cease to account well for the eigenstates behavior at the left of xq. On the 
other hand, a affects mainly the vanishing rate for large x, which affects the 
bound states representation quite severely, but in a rather symmetric way. 
Convergency can be quite substantially improved by tuning the wavefunc- 
tions localization, and this can be achieved by choosing the most appropriate 
s and a, thus precisely the OGQNSB. 

6 Examples of Applications 
6.1 Ar 2 

We compare the OGQNSB and the QNSB for the calculation of the vibra- 
tional spectrum of the Argon dimer, for which a reliable ab-initio molecular 
potential is provided^ in terms of a set of 47 points in the range x = 0.25 to 
20 A. Patkowski et a/™ propose an analytic expression fitting the ab-initio 
points rather accurately. We fit the ab-initio data instead to the expansion 



16 



of Eq. ([T]), up to degree iV max = 8. The resulting best-fit coefficients are 
reported in Table [H Since the repulsive small-a; region does not affect the 
bound states significantly anyway, we privilege the convergence inside the 
binding well region, with a weighted fit.^ Despite its simplicity, generality, 
and the relatively small number of parameters involved (iV max + f = 9), the 
resulting expansion is quite accurate, throughout the whole energy range cov- 
ered by the 47 ab-initio points. In particular, in the well region the agreement 
is quite good, with a RMS discrepancy 5rms of less than half wavenumber, 
see Table [TJ Moreover, the resulting model potential does not suffer from 
the unphysical small-x divergence to — oo of the fitted function,^ and rather 
tracks the repulsive region within few electronvolts. The well depth (classical 
dissociation energy) is D e = ^^> ax (— ) l di = 99.23cm" 1 . 

We apply the algebraic method and solve the resulting quantum-mechanical 
problem ([2]) for the bound-state eigenvalues, using QNSB and GQNSB of 
different size N s . Table [2] compares the results obtained by finite-differences 
solution of the Schrodinger equation for the analytic potential by Patkowski 
et al., 1(> and by numerical diagonalization of the algebraic Hamiltonian (j2J) 
with the parameters from Table[T]on a large N s = 100 OGQNSB (s = 80.18, 
a = 0.213). This large OGQNSB was chosen to ensure that the results 
are fully converged, and is taken as reference. The excitation energies ob- 
tained using our expansion compare favourably to those obtained by using 
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Patkowski et al. analytic expression,^ and to the experimental J = datap^ 
demonstrating equally good or better agreement. 

Table [3] illustrates the convergency properties of the unoptimized QNSB 
by reporting the eigenvalues obtained by diagonalizing the expanded Hamil- 
tonian (TSJ) on N s = 100, 20 and 15 states respectively. The energy differences 
with respect to the N s = 100 OGQNSB reference are shown, in parentheses, 
when exceeding 10~ 3 cm -1 . Fairly well converged results are obtained even 
for the small N s = 15 QNSB. Notice however that the bound state closest 
to dissociation is unbound for N s = 15 and 20, since it is so extended that a 
rather large QNSB (N s > 42) is needed to obtain it at negative energy. Even 
the very large N s = 100 QNSB does not provide a well-converged result for 
that specific level. 

By diagonalizing the expanded Hamiltonian (j2J) on N s = 20 and N s = 15 
OGQNSB, we obtain the complete spectrum, and with an accuracy A — A 
of 6 ■ 10~ 5 and 0.013 cm -1 respectively. The accuracy of all bound levels 
but the last one is basically the same as for the corresponding QNSB, but 
the complete discrete spectrum is obtained, including the highest level. The 
accuracy of the N s = 15 OGQNSB is therefore better than that of N s = 100 
QNSB, for Ar 2 . Reducing the basis size below N s = 15, the highest state is 
missing, but the GQNSB can still be tuned to obtain a fair accuracy of all 
other states (A - A < 0.5 cm" 1 for N s > 11). 
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6.2 H 2 

In a previuos work 1 ' we applied the QNSB formalism to the ab-initio 
adiabatic potential 13 14 for the H 2 molecule. We found that an expansion <^ 
up to iV max = 12 fits all 169 available ab-initio points with a deviation Srms — 
5.5 cm -1 . This expansion, whose parameters are reported in Table IV of 
RefP} produces all the 15 vibrational bound states of this molecule. The 
QNSB parameters for this potential are s = 25.56 and a = 0.564. The 
QNSB produces a cm -1 converged spectrum using N s > 28 basis states. 

By minimizing A, Eq. ( 1241) . we generate an OGQNSB of smaller N s . For 
the calculation of A we use the fully converged N s = 200 QNSB results as 
reference, reported in the second column of Table HI A N s = 25 OGQNSB 
with s = 26.36 and a = 2.115 (A = 0.173 cm" 1 ) produces eigenvalues with 
the same cm -1 figures, i.e. the same accuracy of the N s = 28 QNSB: since 
they are identical to the second column of Table HI they are not shown. For 
less strict accuracy requirements, one could reduce the basis size: the last two 
columns of Table H] compare the eigenvalues obtained with N s = 21 QNSB 
and OGQNSB. The H2 potential expansion illustrates the robustness of the 
GQNSB in state-poor situations: here, for the N s = 21 QNSB eigenvalues 
the differences with respect to the fully converged values reach hundreds of 
wavenumbers, with a RMS discrepancy A = 323 cm -1 , while the discrep- 
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ancy of the eigenvalues obtained by diagonalizing on the N s = 21 OGQNSB 
amounts to A = 3.7 cm -1 only. 

7 Conclusions 

The substantial improvement of the variational accuracy of the bound- 
state spectra computed on a OGQNSB w.r.t. the unoptimized QNSB permits 
in practice to make calculations of a given accuracy on a significantly smaller 
basis size. While this improvement is practically irrelevant to the solution of 
the 1-dimensional vibrational problem of diatomics, it is of great importance 
for the application of this method to the calculation of the spectra based on 
the ab-initio multi-dimensional potential surfaces of polyatomic molecules, 
as is currently pursued in quantum chemical research.^^^^ 1 ^ 1 ^ 1 We are 
currently testing the generalization of the expansion (DO) to the polyatomic 
case 
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TABLE 1: Fit quality and parameters for model potential (PQ), iV max = 8, 
to the Ar 2 potential. 



orms 


[Ha] 


0.26 


3rms 


well [cm -1 ] 


0.48 


a 




0.516787 ao 1 


0-2 




1359.70868 /iHa 


a 3 




1136.96625 /xHa 


a.\ 




181.96578 /iHa 


a 5 




43.51541 /iHa 


a 6 




3.77230 /iHa 


a 7 




-0.13914 //Ha 


a 8 




0.00202 //Ha 






7.116 a 
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TABLE 2: Energy differences (in cm 1 ), between consecutive J = 
vibrational levels of Ar2. 



i-i' 


numerical 2. 


OGQNSB 
N s = 100 


experiment^ 


1-0 


25.76 


25.64 


25.69 


2-1 


20.49 


20.43 


20.58 


3-2 


15.44 


15.46 


15.58 


4-3 


10.79 


10.90 


10.91 


5-4 


6.75 


6.92 


6.84 


6-5 


3.56 


3.71 




7-6 


1.36 


1.31 





"Numerical diagonalization of the Patkowski et al.'s potential.^ 
^Ultraviolet laser spectroscopy data by Herman et al 



TABLE 3: Bound-state eigenvalues of the Ar 2 dimer, computed with 
different methods, all based on the ab-initio values^ [in cm -1 ]. 



state 


OGQNSB 


QNSB 


QNSB 


QNSB 




N s = 100 


N s = 100 


N s = 20 


N s = 15 





-84.41 


-84.41 (-) 


-84.41 (-) 


-84.40 (0.005) 


1 


-58.77 


-58.77 (-) 


-58.77 (-) 


-58.75 (0.01) 


2 


-38.34 


-38.34 (-) 


-38.34 (-) 


-38.31 (0.02) 


3 


-22.88 


-22.88 (-) 


-22.88 (-) 


-22.84 (0.04) 


4 


-11.98 


-11.98 (-) 


-11.98 (-) 


-11.93 (0.05) 


5 


-5.06 


-5.06 (-) 


-5.06 (-) 


-5.02 (0.04) 


6 


-1.35 


-1.35 (-) 


-1.35 (-) 


-1.32 (0.02) 


7 


-0.036 


-0.015 (0.02) 


0.022 (0.06) 


0.051 (0.09) 
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TABLE 4: H2 bound-state energies in reduced-size algebraic bases; in 
parentheses, the differences w.r.t. the reference [in cm -1 ]. 



state 


QNSB° 


QNSBl 


OGQNSET 




N s = 200 


N s = 21 


N s = 21 





-36113 


-36113 (-) 


-36113 (-) 


1 


-31948 


-31948 (-) 


-31948 (-) 


2 


-28020 


-28019 (1) 


-28020 (-) 


3 


-24324 


-24322 (2) 


-24324 (-) 


4 


-20856 


-20845 (11) 


-20856 (-) 


5 


-17614 


-17573 (41) 


-17614 (-) 


6 


-14599 


-14486 (113) 


-14598 (1) 


7 


-11815 


-11588 (227) 


-11814 (2) 


8 


-9271 


-8907 (365) 


-9269 (3) 


9 


-6979 


-6486 (493) 


-6975 (4) 


10 


-4955 


-4376 (579) 


-4951 (4) 


11 


-3222 


-2626 (596) 


-3218 (4) 


12 


-1810 


-1281 (530) 


-1806 (5) 


13 


-761 


-389 (372) 


-757 (4) 


14 


-135 


-9 (126) 


-125 (10) 



"Reference fully converged calculation. 

^Eigenvalues obtained with a N s = 21 QNSB. The maximum difference of 596 cm -1 
corresponds to 1.6% of the well depth. 

Eigenvalues obtained with a N s = 21 OGQNSB (s = 23.52, and a = 3.194). The 
maximum difference of 10 cm -1 corresponds to 0.03% of the well depth. 
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Figure 1. QNSB wavefunctions for n = 0, n = 4 and n = 8, compatible 
with a Morse problem characterized by x$ = 1, a = 4/xq, a 2 = 625, in units 
where fj, — 1, h — 1, so that s = 8.34, a = 0.34. 



Figure 2. Variation of the n = 4 GQNSB wavefunction for a 50% increase 
in the parameters s (a), or a (b), solid line, with respect to the QNSB 
starting wavefunction (dashed line) corresponding to s = 8.34, a = 0.34, 
xq — 1, a — 4/xo, like in Fig. [H 



Figure 3. Discrepancies (E n — E^)/a2 of the individual eigenvalues n 
for the potential V(x) = a 2 [(v(x)) 2 + (w(a;)) 4 ]. (a) N s = 30 (OGQNSB 
with s min = 20.01, a min = 0.435); (b) N s = 16 (OGQNSB with s = 14.47, 
a = 0.314), for all bound states. The OGQNSB (diamonds) discrepancies are 
compared to those based on the QNSB (circles), and to the GQNSB based 
on the choice of (s, a) made by Tennyson and Sutcliffe^"^ (triangles). 



Figure 4. s dependence of A Eq. ( |25|) . for V(x) defined in Eq. (|26l) . 
computed with the GQNSB of N s = 30 elements, as a function of s, and for 
a fixed to o" min , to cr min ± 0.005, and to <r min ± 0.01. 
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Figure 2. 
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